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CN Abstract 

We find a new holographic description of strongly coupled defect field theo- 
ries using probe D5 branes. We consider a system where a large number of probe 
branes, which are asymptotically D5 branes, blow up into a D7 brane suspended 
^ in the bulk of anti-de Sitter space. For a particular ratio of charge density to 

I— I external magnetic field, so that the Landau level filling fraction per color is 

^ equal to one, the D7 brane exhibits an incompressible charge-gapped state with 

p ^ one unit of integer quantized Hall conductivity. The detailed configuration as 

well as ungapped, compressible configurations for a range of parameters near 
the gapped one are found by solving the D5 and D7 brane embedding equations 
numerically and the D7 is shown to be preferred over the D5 by comparing their 
^ energies. We then find integer quantum Hall states with higher filling fractions 

as a stack of D5 branes which blow up to multiple D7 branes where each D7 
O brane has filling fraction one. We find indications that the u D7 branes describ- 

ing the filling fraction u state are coincident with a residual SU(z/) symmetry 
when is a divisor of the total number of D5 branes. We examine the issue of 
stability of the larger filling fraction Hall states. We argue that, in the D7 brane 
CS| phase, chiral symmetry restoration could be a first order phase transition. 



1 Introduction and discussion 

The D3-D5 brane system is one of the best studied examples of top-down holography 
using probe branes [I]- [20]. A stack of A^5 coincident and appropriately oriented D5 
branes is placed in the vicinity of a stack of coincident D3 branes. In the limit 
where the number of D3 branes N is taken large and the string coupling gs small 
so that A = 4:7rgsN is held constant, the D3 branes are replaced by the AdS^ x 
background spacetime and the D5-branes are a probe embedded in this background. 
The probe limit requires that « N. The string theory is weakly coupled in 
the limit where A is large. In this limit, the mathematical problem of understanding 
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the probe branes reduces to understanding the embedding of their worldvolume in 
AdS^ X S\ 

The field theory dual of the large limit of string theory on AdS^ x is the planar 
limit of A/" = 4 supersymmetric Yang Mills theory with 't Hooft coupling A. The role of 
the probe branes is to introduce fundamental representation degrees of freedom. In the 
D3-D5 system, these degrees of freedom constitute a supersymmetric hypermultiplet 
which occupies a 2+1-dimensional defect embedded in 3+1-dimensional spacetime. 
This system has a solution with superconformal symmetry which preserves half of 
the supersymmetries of the A/" = 4 theory and the conformal symmetry of the defect. 
The probe brane worldvolume is the space AdS^ x which has bosonic symmetries 
5^0(3, 2) X S0{?>) X S0{?>). When there are coincident D5 branes, there is a U{N^) 
fiavor symmetry in the field theory dual. 

When temperature, an external magnetic field and charge density are introduced, 
the D3-D5 system has a rich phase diagram [13j. For example, at zero temperature 
and density, the presence of a constant external magnetic field results in spontaneous 
breaking of an SO (3) chiral symmetry of the dual field theory. This is an example of 
magnetic catalysis of chiral symmetry breaking, where the tendency of interactions 
to break chiral symmetry is greatly enhanced in the presence of an external magnetic 
field [21j-[32j. Due to chiral symmetry breaking, the field theory spectrum gets a 
charge gap. In the holographic theory, this is seen as the probe brane taking up a 
"Minkowski embedding", where it still approaches AdS^ near the boundary of AdS^^ 
but it is modified in the bulk and it caps off before it reaches the Poincare horizon. 
The charge gap is seen in the fact that the D3-D5 strings, whose low energy modes 
are the U(l) charged fundamental representation matter degrees of freedom, must 
be suspended between the D5 worldvolume and the Poincare horizon. When the 
worldvolume does not reach the horizon, the D3-D5 strings must have a minimum 
length and therefore they have an energy gap. 

When a charge density is introduced in addition to the magnetic field, there is a low 
density phase which still exhibits dynamical chiral symmetry breaking. In all cases, 
when the charge density is nonzero, the D5 brane embedding must be a "black hole 
embedding", rather than Minkowski embeddingj^ It must stretch from the boundary 
of AdS^ to the Poincare horizon. This is due to the fact that it carries worldvolume 
electric fiux. It is not possible for the worldvolume to pinch off smoothly unless 
there is a sink to absorb the electric fiux. Fundamental strings suspended between 
the worldvolume and the horizon could provide such a sink. However, fundamental 
strings with the necessary configuration always have a larger tension than the D5- 
brane, and they would pull the D5 brane to the horizon. Indeed, the numerical 
solutions for the worldvolume geometry do exhibit a spike which resembles such a 
fundamental string funnel. A black hole embedding is thus a compressible state. It 
does not have a charge gap. It is somewhat puzzling that the D5 brane states do not 

^The term "black hole embedding" which normally used to describe a probe brane with worldvol- 
ume which reaches an AdS^ black hole horizon as well as intersecting the boundary of AdSr^^ Here 
we shall use this term, even though we shall deal solely with the zero temperature limit where there 
is no black hole. We use it to refer to probe brane whose worldvolume reaches the Poincare horizon 
of the Poincare coordinate patch of AdS^. 
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have a charge gap and are compressible for any non-zero value of the charge density. 
In this paper, we examine the question whether the U(l) charged D5 brane can have 
an incompressible state at all. We shall find that, for an interesting range of U(l) 
charge, the answer is yes and we will construct the corresponding solutions. 

In rough outline, our solution of this problem is as follows. We consider A^5 D5 
branes which enter AdS^ at its boundary. At the boundary the D5 brane worldvolume 
wraps a 2-sphere S'^ embedded in and has an asymptotically AdS^ component 
embedded in AdS^. It also has worldvolume gauge fields which are needed so that the 
dual field theory has a constant external magnetic field and a constant U{1) C U {N^) 
charge density. For simplicity, we will consider only the case where the temperature 
is zero, although we anticipate that most of our results could easily be generalized 
to finite temperature. As the AdS radius decreases, the D5 branes blow up via the 
Myers effect [33j to a D7 brane whose worldvolume wraps a second two 2-sphere 
in S^. The second 2-sphere carries a Dirac monopole bundle of the D7 worldvolume 
gauge field with A^5 units of monopole charge. The same coupling of the D7 brane to 
the Ramond-Ramond 4-form fields which implements the Myers effect also generates 
a topological term ^ F A F A C^"^^ in the D7 brane worldvolume gauge theory. This 
term allows the magnetic fiux to carry electric charge. In this way, the electric charge 
dissolves into the D7 brane. For a certain value of the charge density, all of the charge 
can be dissolved and the D7 brane is allowed to cap off at a finite AdS radius, that is, 
to have a Minkowski embedding. The result is an incompressible state with a charge 
gap. The filling fraction (the dimensionless ratio of charge density to magnetic field) 
at which this occurs is precisely that of a. u = 1 integer quantum Hall state. If the 
filling fraction u is close to, but not equal to one, the D7 brane still exists, but it has 
a different character as it a has residual charge and it must therefore have a black 
hole embedding which reaches the Poincare horizon. We are able to confirm this by 
numerical solution of the embedding equations for the D7 brane where we find the 
Minkowski embedding at = 1 and black hole embeddings for the D5 brane with any 
nonzero charge density and the D7 brane when u 1. We also find that these D7 
brane solutions have lower energy than D5 brane solutions with the same parameters, 
and are therefore favored, at least in an interesting parameter region including u = 1 
and perhaps a much larger region which is yet to be explored. Our analysis is most 
reliable when A^5 ^ a/A which is large. Of course, our holographic computations are 
reliable when \/A is large. 

Then, we also find higher integer quantized Hall systems by considering, for = 2, 
for example, A^5 D5 branes where a subset of them, say N^a D5 branes, blow up 
into one D7 brane and the remaining N^i^ = — N^a D5 branes blow up into a 
second D7 brane and where both of the latter branes are in the u = 1 charge-gapped 
incompressible state. For = 2, we show that this state is preferred over both the 
un-gapped u = 2 D7 brane solution and the D5 brane with the same parameters. We 
also find that the symmetric solution where N^a = A^56 = A^s/S seems to be preferred. 
This state has a residual SU (2) symmetry, and to have this symmetry in the strictest 
sense, A5 must be an even number. We speculate that there can be a large number 
of stable integer quantum Hall states up to the maximum filling fraction which can 
be achieved in this way, u = although many are outside of the domain where our 
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Figure 1: The Landau level spectrum of an S0(3) doublet of massless non-interacting fermions 
in a magnetic field is depicted. The energy levels are the horizontal lines and the vertical axis is the 
single-particle energy. Only fermionic levels are shown here. We have assumed that A^5 = 1 and we 
have suppresed the color degeneracy. The two columns represent the two isospin states. Fermions 
are depicted by the dots. In the neutral ground state, negative energy levels are filled and positive 
energy levels are empty. The zero energy states must be half-filled. In the absence of interactions, all 
possible fillings of the zero energy states have the same energy. This gives rise to a large degeneracy 
of the many-body ground state. 

analysis using probe branes is reliable. In order for a higher integer filling fraction to 
be SU{v) symmetric, v must be a divisor of N^. 

1.1 Weak coupling 

It is interesting to examine what one would expect in the weak coupling limit. Our 
purpose here is to understand the origin of Hall states and the role of the scalar field 
which is necessarily present in this supersymmetric theory. If the coupling constant 
is set to zero, the action of the defect theory is for a A^5 D5 branes is known to be ^5] 

Sb = j d''x{\D,q\^ -i^YD^^], 

where g is a complex scalar field which is a spinor representation of /SO (3), the 
isometry group of the worldvolume aS^, and is a complex fermion which is in a 
spinor representation of S'0(3), the isometry group of S'^ as well as a spinor of the 
5^0(2, 1) Lorentz group. The 2-sphere S'^ is the one which was originally wrapped by 
the D5 brane and which the D7 brane worldvolume also wraps. The 2-sphere S'^ is 
the one which is formed as a fuzzy sphere of D5 branes and which is wrapped by the 
D7 worldvolume and contains A^5 units of D7 worldvolume magnetic flux. Both q and 
^ transform in the fundamental representation of the bulk SU(N) gauge group and a 
flavor SU(A^5) symmetry. They have a color index which runs from 1 to N and flavor 
index which runs from 1 to A^5. 

It is easy to see that, in the presence of a magnetic field, this system will break the 
aS0(3) symmetry when there is even an infinitesimally weak interaction. In an external 
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U{1) magnetic field, the single particle energy levels of the fermions are the Landau 
levels En = \f2Bn and the scalar field are uj^ = \/ i^n + X)B with n = 0,1,2,... in 
both cases. The levels are infinitely degenerate. The density, i.e. the number of states 
per unit volume, of a single level is The degeneracy of each fermionic Landau 
level is then 2N^N^ where the factor of 2 is due to the fact that the fermion is an 

5^0 (3) spinor, A^5 and come from the numbers of fiavors and colors, respectively. 

The bosons have a gap and are not excited at energies less than (jJq = \fB. The 
fermions, on the other hand, have a zero energy mode, which will dominate the 
low energy properties of the system. Moreover, in the neutral ground state of the 
system, depicted for the case of one fiavor, A^5 = 1 (and two S'0(3) isospin states), in 
figure 1, all negative energy fermion states must be occupied, positive energy states 
must be empty and the zero modes should be half-filled. The latter requirement is 
due to a charge conjugation (or particle-hole) symmetry of the spectrum of the Dirac 
Hamiltonian in a magnetic field. 

An arbitrarily weak interaction will have significant infiuence on the zero modes. 
The perturbative shift of the energy of states such as in figure [T] due to a weak 
repulsive interaction like the Coulomb interaction is known to be minimal for the 
states depicted in figure [2} where all of the zero mode fermions take up one of the 
S'0(3) spinor indices [26j. This is similar to, and has the same source as the Hund 
rule for populating electron energy levels of atoms. In rough terms, an interaction 
like the Coulomb interaction, which is repulsive at all scales, would be minimized 
by keeping the fermions as far apart as possible. Since their wavefunction is totally 
antisymmetric, this is achieved by making it as antisymmetric in particle positions as 
possible, and therefore as symmetric as possible in the other indices. Symmetry in 
the fiavor indices is maximal if all states have the same index. This is the symmetry 
breaking state depicted in figure [2] This phenomenon is known in the condensed 
matter physics literature as quantum Hall ferromagnetism and it is the explanation 
of an experimentally observed feature of some two dimensional electron gases [34]-[37]. 

Now, assume that interactions have generated a small gap. First, we observe that 
this breaks the 5^0(3) chiral symmetry and it makes the neutral state a gapped state, 
where otherwise, if the symmetry was not broken, it would have been ungapped. 
What will the system look like if we now vary the density? If we begin with the 
neutral state of figure [2] and add fermions we immediately have an ungapped state 
with the second set of would-be zero modes partially filled. Then, if we continue to 
add fermions, we will eventually fill all of the remaining would-be zero modes which 
are now at a small positive energy. Then the completely filled state, depicted in figure 
[3} also has a gap. This state is an integer quantum Hall stat^ with filling fraction 

^Note that, in this translational invariant quantum field theory, even though we have gapped 
states for certain discrete values of the charge density, the usual signature of the quantized Hall effect, 
the Hall plateaus in the graph of Hall conductivity versus filling fraction will not occur. Of course, 
localized states, which are absent here, would be needed in order for integer quantum Hall plateaus 
to form and for this system to have a conventional quantum Hall effect. Here, we nevertheless refer 
to the charge gapped state obtained by tuning the charge density to the appropriate value as a 
"quantum Hall state". It is distinguished from other states by virtue of being incompressible and by 
the existence of a charge gap. 
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Figure 2: Magnetic catalysis: with an arbitrarily weak repulsive interaction (such as the Coulomb 
interaction), the preferred state is one where all of the zero modes with one label are occupied to 
obtain the u = state. This is the analog of the Hund rule for electrons filling atomic orbitals. 
It can be shown to minimize the exchange energy. This is a gapped state which, since the density 
vanishes, has vanishing Hall conductivity. 



1, where 1^ 



The color degeneracy {SU (N) of A/" = 4 Yang-Mills theory) is explicit in our definition 
of filling fraction. 

If we continue to increase the density, the next available energy states are where 
we excite bosons whose gap is y/B. Then there can be no further Hall states as all 
further charge can be accomodated by Bosons which do not have a Pauli principle. 
If we decreased, rather than increased the charge density, particle-hole symmetry 
tells us that we would also find an incompressible state at = — 1. Thus, at weak 
coupling, when there is only one flavor of fermions, and because of the presence of 
charged scalar fields, there are only three incompressible states with = 1, 0, — 1 and 
only two Hall states with = 1, — 1. This changes if there are multiple flavors. 

For multiple D-branes with A^5 flavors of fermions, there are similar arguments 
to quantum Hall ferromagnetism for a complete resolution of the degeneracy of the 
zero modes so that filling an SO (3) state for each flavor (with fermins or holes) could 
potentially yield 2N^ quantum Hall states. However, the number of these states which 
actually become Hall states depends on how many fit under the threshold for creating 
a scalar particle. At very weak coupling the number is A^5 levels with filling fractions 
±iy = 1,2,...,A^5 and Hall conductivities ±o-xy = ^5^5 An example with 

A^5 = 3 is depicted in figure |4} 



^Note that this filling fraction is defined with a factor of the number of colors N in the denomi- 
nator. The first intger filling fraction u = 1 refers to a Landau level which is completely filled with 
baryons, where the baryon is a composite of N fundamental representation quarks and has U(l) 
charge N. This color degeneracy is unavoidable here and it is natural if the theory is confining, as 
one would expect when the chiral symmetry is spontaneously broken, as it is here. 
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Figure 3: Beginning with the state depicted in figure [2] and increasing the chemical potential 
and charge density eventually results in the u = 1 state that is depicted here. The energy needed 
to excite the first bosonic states is depicted by the dashed line. Further increasing the charge will 
simply increase the number of Bosons excited. Thus, the u = 1 state is the only quantum Hall state 
in the spectrum. 
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Figure 4: For = 3, the possible Hah states have filling fractions u = 1,2,3. The state with 
= 2 is depicted. There is an energy gap between the fermions in this state and those which must 
be added to continue filling the remaining zero modes. 
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In effective field theory, the HaU conductivity is encoded in a Chern-Simons term. 
If, in one of the gapped states, we introduce a background U(l) gauge field and then 
we integrate out the matter fields, the action has a Chern Simons term 

Ses. = ^j£x e^^^A^d^Ax + ... 

When there is a charge gap, the Chern-Simons term does not renormalize beyond one 
loop [38] [39] [40] . It is also known that when there is no gap, it can be corrected at 
higher loops [39]. At one loop, the scalar field does not contribute and the fermion 
gives the result above [41][42j. 

The no-renormalization theorem suggests that the filling fraction and Hall con- 
ductivity of a gapped state could be exact and identical to its value in the weak 
coupling limit. It is interesting to check this idea at strong coupling. In this paper, 
we indeed find, first of all, the Hall state with the lowest integer Hall conductivity at 
filling fraction u = 1^ and then others with = 2, 3, . . . the total number limited by 
dynamical issues, but bounded from above to be less than or equal to N^. 

One might speculate that these states are strong coupling extrapolations of the 
ones that are seen in the weak coupling limit. It is interesting that the upper limit on 
the strong coupling side comes from the D brane geometry. We shall also see that the 
maximum number of Hall states is A^5, coming from A^5 D5 branes blowing up into 
A^5 D7 branes. However, it is very possible that only a much smaller number of these 
states are stable. The geometric description of the blow-up as a classical sphere with 
magnetic fiux, which our analysis relies on, is also not accurate unless the number of 
D5-branes which blow up to a D7 is large. 

1.2 Dr 

The D7 brane configuration that we shall use in this paper is closely related to that 
used in a series of papers [^-[48] where one examines a D3-D7 system where the D7 
brane worldvolume wraps S'^ x S'^ C and the D7 brane is stabilized by putting 
U(l) worldvolume magnetic fiux on S'^. This configuration is called D7' and it has 
been used to construct a model of the Hall effect by a mechanism which has some of 
the ingredients of the one that is used here. Both D7 and D7' probe branes have also 
been used to model some phenomena in planar condensed matter systems [50]- [55]. 

An important feature of a Hall state is that it is incompressible, that is, it has an 
energy gap for charged excitations. For a probe brane to be incompressible, it should 
have a Minkowski embedding. However, to be a Hall system, it must also have electric 
charge and at first sight these are incompatible, as we have discussed above, when 
a probe brane carries the worldvolume gauge field configurations necessary for its 
field theory dual to have electric charge, it should have a black hole embedding. A 
way around this occurs in the D7' system exactly as it does in the system that we are 
interested in in this paper and it was already used to construct Hall states in reference 
[13]. In the presence of external magnetic field, the Wess-Zumino terms in the D7 
brane action allow electric charge to be dissolved into the brane worldvolume. This 
could occur for discrete values of the ratio of charge density to magnetic field, the 
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filling fraction. Then, in these special cases, the solution for the D7 brane worldvolume 
is a Minkowski embedding with a charge gap. These were identified as quantized Hall 
states of D7^ Due to a different choise of boundary conditions, the filling fractions 
are not integers in that case. With the types of embedding chosen for some 
internal magnetic fiux on the D7 brane worldvolume was required for stability of the 
system. The fiux was also quantized and that gave a quantization condition for Hall 
conductivity (which did not match either the integer or known fractional quantum 
Hall effects). 

The difference between D7' and what we do in the present paper is that we choose 
a different boundary condition at the boundary of AdS^. The boundary conditions 
at the boundary of AdS^ must be chosen so that the numerator of the last term 
in the differential equation which describes the behavior of the 5-sphere latitude at 



which the D7 brane is embedded, equation (35), vanishes there (at large r in (|35|)). 
The numerator of that last term is the partial derivative of the expression in the 
denominator whose large r behavior is the 'potential' 

Viip) = 4 sin^ V^(/^ + 4 cos^ ip) . 

We thus need the asymptotic value of ip{r) to approach an extremum of this po- 
tential. Moreover, if we want the fiuctuations about the asymptotic value to have 
two normalizable modes, it is necessary that this extremum is a maximum. Simple 
calculus tells us that there are either one or two maxima: if < |, there are two 

maxima at = | arccos ^1 — 2p — 1 ^ and ip = ^. On the other hand, if 

> ^, there is only one maximum at = ^. In reference [43], they used the first 



maximum. In the present paper, we will consider of the second maximum at = |. 
This as a place where cos^ = and the D7 brane world- volume degenerates to a six 
dimensional hypersurface where one of the worldvolume spheres, 5"^, has shrunk to 
zero volume. This was a sphere carrying magnetic fiux. When it shrinks to a point 
it leaves a singularity in the form of a point Dirac monopole. This singularity should 
be interpreted as A^5 = ^/ D5 branes ending at the D7 brane. In fact, we interprete 
it as a region where a number A^5 of D5 branes takes on a non-Abelian configuration 
and blows up to form the D7 brane. 



1.3 Outline 

In section 2 we shall review the problem of determining the geometry of a D5 brane 
embedded in AdS^ x when the field theory which is its holographic dual has an 
external magnetic field and charge density. In section 3 we discuss the embedding of 
the D7 brane with the same external magnetic field and charge density. In section 4 
we discuss our numerical solutions for the D5 and D7 brane configurations as well as 
their energies. Section 5 contains some conclusions. 
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2 D3-D5 



We shall study the D3-D5 system in the probe limit where the number of D5 branes 
A^5 is much smaller than A^, the number of D3 branes. The material of this Section 
is already well known and can be found elsewhere in the literature. We include it 
here for the convenience of the reader, and to fix our conventions and notation which 
will be needed later on when we compare the D5 brane with the D7 brane. On 
fiat space, the D3 and D5 branes are oriented as in Table 1. They overlap in 2+1 
dimensions. In principle, with multiple D5 branes, their coordinates are matrices and 
the worldvolume gauge fields have U {N^) gauge group. We will begin with an analysis 
which assumes that the non-Abelian structure plays no role. The actions for each of 
the individual D5 branes are identical and it sufiices to study one of them and to 
multiply the total action by the number of branes A^5. 
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Table 1 : D3, D5 and D7 orientation 

We shall be interested in the limit of the string theory which coincides with the 
planar limit of the gauge theory and, after the planar limit is taken, the large A 
strong coupling limit. In this limit, the string theory is classical, and the problem 
of including a D5 brane in the AdS^ x geometry reduces to that of finding an 
extremum of the Dirac-Born-Infeld and Wess-Zumino actions, 



S = -N, 

9s 



V-det(£/ + 27rtt'J') + 2t:q'C^^^ A J" 



(2) 



where Qs is the closed string coupling constant, which is related to the A/" = 4 Yang- 
Mills coupling by ^tiQs = Qym-) ^he coordinates of the D5 brane worldvolume, 
gab{o') is the induced metric of the D5 brane, C^^^ is the 4- form of the AdS^ x 
background, is the worldvolume gauge field and 



(3) 



is the D5 brane tension. Note the overaU factor of A^5 in equation The Wess- 
Zumino action wiU not contribute to the D5 brane equations of motion for the types 
of embeddings that we will discuss here. 

We shall work with coordinates where the metric of AdS^ x is 



ds^ = VA, 



a 



'^{-df + dx^ + dy'^ + dz^) + ^+ 



(4) 
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Here, {t^x^y^ z^r) are coordinates of the Poincare patch of AdS^. The boundary of 
AdS^ is located at r ^ oc and the Poincare horizon at r ^ 0. The coordinates of 
are a fibration of the 5-sphere by two 2-spheres over the interval ij) G [0, 7r/2]. 

The dynamical variables are the ten functions of six worldvolume coordinates 
which embed the D5 brane in AdS^ x S^^ as well as the six worldvolume gauge fields. 
We shall take the static gauge where the D5 brane has coordinates (t, x, ^, r, 0, 0) 
and, due to symmetry, the remaining two coordinates z{r) and ip{r) depend only 
on the AdS radius r. The equation of motion for z is satisfied by a constant. The 
worldvolume metric is 



ds^ = VXi 



a 



r'^{-dt^ + dx^ + dy^) + 




(5) 



If rj^^{r) ^ at r ^ oc, in the large r regime, the worldvolume is a product of 
AdS^ and . We shall also take the worldvolume gauge fields with field strength 



a 



d 
dr 



a{r)dr A dt + bdx A dy 



The constant external magnetic field is 

and the r-dependent temporal component of the gauge field is 

A{r) = ^ a{r). 



(6) 



(7) 



(8) 



With the above ansatz for the worldvolume metric and gauge fields, the Born-Infeld 
action becomes 



^5 = -ATs A^5 dr2 sin^ t/^V^^ + Hy^l + 
where, using ([s]), 



dip 
dr 



da 
dr 



(9) 



2VXN 
(27r)3 



2+1- 



(10) 



The factor of — in (10) is the coefficient of the Dirac-Born-Infeld action (2), T5 is 

9s ] _j ] [ 

rom the 



given in ([3|) and we recall that gs = A/(47rA^). The factor (a/A<^0^ comes 
overall factor in the worldvolume metric in equation The factor of (27r) is from 
the integral over the worldvolume two-sphere. The integral over (x, t) produces 
the volume factor V2+1. 



^It is half of the volume of the 2-sphere. The other factor of 2 is still in the action in front of 
sin^ t/j. This notation is designed to match with the D7 brane, which we shall study in the next 
section, and to coincide with notation in reference [43j . 
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Now, we must solve the equations of motion for the functions i/j^r) and a(r) which 
result from the Lagrangian ^ and the variational principle. We will use the boundary 
condition 

TT 

lim i){r) = -, (11) 

which is compatible with the equation of motion for ip. Since the Lagrangian depends 
only on the derivative of a(r) and not on the variable a(r) itself, a(r) is a cyclic variable 
and can be eliminated using its equation of motion. 



d SS, _^ 2sin2^V^(62 + r4)|a 



_ Q _^ ^ V ^ ^ dr ^ (12) 

where is a constant of integration. We can solve for ^^a. 



— g = ] (13) 

dr ^4sinV(&^ + ^^) + ?5 

We now must solve the remaining problem of determining ip^r) where we wish to 
extremize the Dirac-Born-Infeld action for a fixed value of the integration constant 
^5 which will be proportional to the total electric charge. To this end, we use the 
Legendre transformation 



to eliminate ^a. We obtain the Routhian 

ar 



7^5 = -M,N, J rfrY/4sinV(&' + r4) + 95^1+ , (14) 

which must now be used to find an equation for '^(r). Applying the Euler-Lagrange 
equation to the Routhian ([T4|) leads to 



1 + (rl^) dr 



8r sin ih 

1 + 



Asin'^ xp{b'^ + r*) + ql 



8 sin^ ip cos -0(6^ + r"^) 
4sinV(^^ + ?'^) + ?5 



0. 



Remembering that the total charge density is defined by 

_ 1 SS5 _ 1 2n SS5 _ 1 27r _ 2NN5 

we see that is related to the total charge density p by 
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Without loss of generality, we can rescale r in the equation of motion so that 6=1. 
Then, in that equation is replaced by ^5/6, which, for future reference we write as 



b 



T' 



(16) 



where the parameter u defined in ([T]) is the filling fraction of a free particle Landau 
level (the density of a Landau level is 5/27r), normalized by a factor of the number of 
colors of the bifundamental charged fields. If of the fundamental "quark' fields must 
bind to form a baryon, which has U(l) charge A/', then z/ would be the filling fraction 
of a single Landau level of noninter acting baryons. The factor in the denominator of 



(16) is related to the total number of D5 branes, 



(17) 



This will be an important parameter in the next section when we consider the D7 
branes. Here we note that we are interested in the parameter regime where / is of 



order one when A is large. This means that the number of D5 branes, 
large. Finally, we present the equation of motion as 



- Va 

27r 



/is 



+ r—ip 
ar 



1 



8r* sin*^ i/jf"^ 



4sinV/^(l + r^) + (ttz^)^ 



4smV/^(l + r4) + (7ri/)2 



0. 



(18) 



The Routhian from which this equation is derived is 



f \Vx 



dr \ Usm'^xpP{l + r*) + {iruy 




(19) 



We can now see that the boundary condition that t/^ ^ | at r ^ (X) is compatible 
with the equation of motion. In fact, ^0 = | is a solution of the equation of motion for 
all values of r. However, this constant solution is known to be unstable if the filling 
fraction is small enough [12j. The critical value of the filling fraction is 



= 2y7//7r - 1.68/. 



(20) 



When \u\ < v^) the stable solution of the equation of motion has r-dependent ij). 



When |z/| > z/c, constant ifj 
ip{r) in the large r regime. 



dr 



is the stable solution. Looking at the equation for 



d_ 
dr 



5^ + 25^ = 0, 



we see that it must have the asymptotic behavior 

TT Cl 
- + — 

2 r r 



'0(r) 



C2 
2 



(21) 



(22) 
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The constant Ci has the holographic interpretation of being proportional to the bare 
mass of the fundamental representation field^ The constant C2 is proportional to 
the chiral condensate. 

One of these constants is to be specified in the boundary conditions which define 
the theory. The other constant is then determined by solving the equations of motion. 
In this Paper, we will restrict our attention to the massless theory, that is, we will 
set Ci = and and, consequently, C2 will be determined by the dynamics. 

Then, if it happens that C2 = 0, because ip must satisfy the second order differential 



equation (18), it must be a constant , ^0 = ^ for all values of r. This is the solution 



2 

with vanishing chiral condensate which is stable when the charge density is big enough, 

If C2 7^ 0, the chiral condensate is nonzero and ip is r-dependent. Since the 
bare mass is zero, the chiral condensate being nonzero is interpreted as spontaneous 
breaking of chiral symmetry. This phase is stable when \u\ < v^- 

The phase transition which occurs when = is thus interpreted as a chiral 
symmetry breaking/restoring phase transition which is driven by decreasing/increasing 
the quark density. This phase transition was found in reference [12]. One very inter- 
esting feature is that it has Berezinsky-Kosterlitz-Thouless scaling and it is thus one 
of the few holographic phase transitions that does not exhibit mean field scaling. 

We also note that, when 7^ 0, even when chiral symmetry is broken and the 
quarks have dynamically generated masses, the finite density system that we find is 
compressible. As we have already discussed, the D5 brane must necessarily reach the 



Poincare horizon. One way to understand this from the Routhian (19) is to note that, 
if the brane were to cap off at a finite radius ro, the derivative dip/dr would diverge 
there. This would result in a very large energy unless the first square root factor in 
the Routhian compensated it by going to zero. This cannot happen when 7^ 0, 
since that factor is bounded from below by \7Tiy\. In the neutral case, when z/ = 0, it 
is allowed, t/^ ^ as r ^ Tq and the square root factor would indeed vanish in that 
special case and it can thereby keep the energy finite. 

In the next section, we shall derive the equations of motion for the D7 brane 



with the same charge density, magnetic field and asymptotic boundary condition (22) 
where we will also set Ci = 0. 



^ The constant ci can indeed be seen to be proportional to the asymptotic at large r physical 
separation of the D5 and D3 branes. When this separation is non-zero, the open strings connecting 
the D3 and D5 branes have a minimum energy, or a mass gap. This is interpreted as a bare mass 
for the hypermultiplet fields of the field theory dual. Note that such a mass breaks and SO (3) 
symmetry expHcitly, by separating the D3 and D5 branes in one of the transverse directions. Since 
the SO (3) rotates the transverse directions into each other, specifying the direction of separation 
in those directions breaks the symmetry. The field theory operator dual to fluctuations of this 
separation is a parity symmetric, S0{3) breaking mass term for the fermions in the hypermultiplet. 



The exponents 1 and 2 of ^ and ^ in (22) correspond to conformal dimensions of the mass and 
mass operator which, since they are protected by supersymmetry in the D3-D5 system, coincide 
with their tree level, engineering dimensions. 
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3 D3-D7 



We shall study the D3-D7 system in the same probe limit in which we considered 
the D5 branes in the previous section. On flat space, the D3, D5 and D7 branes are 
oriented as in Table 1. In this case, there will be a single D7 brane whose worldsheet 
will wrap both of the spheres S'^ and S'^. The worldsheet gauge flelds have A^5 units 
of magnetic flux on S'^. The Dirac-Born-Infeld and Wess-Zumino actions for the D7 
brane are 



(23) 



where 



(27r)7a 



7^/4 • 



(24) 



is the D7 brane tension. 

We shall work with coordinates where the metric of AdS^ x are given in ^ and 
the Ramond-Ramond 4-form of the IIB supergravity background in this coordinate 
system is 



r^dt A dx A dy A dz -\ cos 6 Adcj) Ad cos 9 A d(j) 



(25) 



Here, d^c{ip) = 8 sin^ t/^ cos^ t/^ which has the indefinite integral 

1 TT 

c(^) = iIj- - sin 4^ - -. 



(26) 



Here we have chosen the integration constant so that c(7r/2) = 0, i.e. it vanishes at 
the asymptotic value of the angle ip. The choice of this constant is a string theory 
gauge choice and our results will not depend on its speciflc value. 

In the "static gauge" the D7-brane has coordinates (t, x, r, ^, 0, 0, 0). The equa- 
tion of motion for z is satisfled by a constant. The D7 brane worldvolume metric 
is 



ds^ = 



a 



r'^{-df + dx^ + dy'^) + 



1 + r 



+ sin^?/;((i(9^ + sin^ 

The worldvolume gauge flelds are 

d 



') + cos^ilj{d9^ + sin 



/ 




2i{olT = ^f\oL ( ^a{r)dr Adt + bdx Ady + —d cos 9 A 
dr 2 



(27) 



(28) 



he parameter / is the flux of the worldvolume gauge flelds deflned in equation 
(|17). It corresponds to Dirac monopoles on 5"^. The external magnetic fleld is 
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identical to that on the D5 brane. The temporal component of the gauge field will 
be determined so that the total charge density is p. Substituting the ansatz into the 
Dirac-Born-Infeld action, we obtain 



^7 = -M, 



dr 



2 sin^ ^^J{p + 4 cos4 ^) {h^ + W 1 + (r 



dip 
dr 



da 
dr 



, (29) 



where, using (24), 



A/V 



2AA^ 
(2^ 



V2+1. 



(30) 



We shall use the same boundary condition for ip^r) that we imposed in the case of 
the D5 brane in (22). This boundary condition is compatible with the equation of 
motion for ip^r) which we shall derive below. 

Now, in a completely analogous way to the D5 brane case, we must solve the 
equations of motion for the functions ip{r) and a(r) which result from the action (29) 
and the variational principle. Again, a(r) is a cyclic variable and can be eliminated 
using its equation of motion. 



d 5S7 
dr 5-fa{r) 



^ 



2 sin^ ^^{P + 4cos4 7/;)(62 + r'^)£a 



- 2bc = gr. 



(31) 



where qj is a constant of integration. It is proportional to the total charge density in 
the field theory dual 



1 



P = 



5Sj 



1 2n 



(32) 



We can solve for ^a. 



d 



-a = 



dr ^4 sin^ t/;(/2 + 4 cos^ ^) (52 + r^) + (g^ + 26c)2 ' 
The Legendre transformation yields the Routhian 



(33) 



re 

7^7 = -^f^ 

Jo 



drA/4sin^^(/2 + 4cos4^)(52 + r4) + (g7 + 26c)2Wl+ r 



dip 
dr 



(34) 
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which must now be used to find an equation of motion for ip^r). Applying the Euler- 

ieads to 

8r^ sinV(/^ + 4cosV) 



Lagrange equation to the Routhian (34) leads to 



1 + irf) 



+ r— 
^ dr 



1 + 



4sinV(/^ + 4cos4^/;)(l + r^) + {n{p - 1) + 2^ - | sin 4^)2 
8sin^^cos?/^/2(l + r*^) + 4sin^ 2ip cos 2iljr'^ + Asin^ 2'4){'k{p - 1) + 2il}) 



4smV(/^ + 4cos4^)(l + r4) + {'k{p -l) + 2^-\ sin4?/;)2 



(35) 



We have rescaled r so that 6 = 1. Here 

q Ati^ a 1 



6 VAA^ 27rS 
where v is defined in equation ([T]). The Routhian is 



^7 = 



(36) 



where 



Lt = A/4sinV(/^ + 4cos4 7/;)(l + r^) + (7r(zy - 1) + 2^^ - - sin47/;)2W 1 + ( r 



(37) 



Now, we note the similarity of the Routhian (37) for the D7 brane to that of the D5 



brane (|19|). At the asymptotic boundary, r ^ cx), the Lagrangians become identical. 

Af7f=Af5N,. (38) 



To see this, we note the identity 



We can now see that, from the form of the Routhian, the D7 brane can have a 



Minkowski embedding when the first square root in L7 in (37) can go to zero in order 



to compensate the fact that at a smooth cap-oflF, dip/dr must diverge and make the 
second factor large. The first square root contains a sum of squares which can go to 
zero only when ip goes to zero and only when v = 1^ the charge density coinciding 
with the first integer Hall filling fraction. 

We note that the r-independent chirally symmetric solution, = | is also a 



solution of the D7 brane equation of motion (35). We would expect that this solution 



is stable if the charge density v is large enough. In fact, by linearizing the equation 
about the constant solution, we can see that it is stable to small fiuctuations in 
precisely the same region that the D5 brane is stable, for \v\ > v^- When \v\ < v^^ 
like the D5 brane, the D7 brane must have an r-dependent solution which breaks 
chiral symmetry. Our numerical results will indicate that there are r-dependend D7 
solutions even in the region where the constant solution is stable to small fiuctuations, 
and these solutions have lower energies than the constant solution. This means that 
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the chiral symmetry breaking solution is actually the prefered one there, and that 
there is an energy barrier between the two. 

Now, we can compare the mathematical problems of finding the D5 and D7 brane 
embedding with the same values of external parameters, S, v and Ci, the latter being 
set to zero. 



4 Numerical Solutions 



We must solve the equations (18) and (35) with the boundary condition (22) where 



we set Ci = 0. For each differential equation we start by establishing a high order 
power series solution, valid at large r, which obeys the boundary condition and has C2, 
/ and V as parameters. With the aid of these series solutions we then generate initial 
data for i\){r) and r^'0(r), at a large value of r, which we use as input to a numerical 
integration routine in Mathematica that allows us to integrate the prospective solution 
to smaller values of r. In this process we use a shooting technique which, for a given 
set of values of (z/, /) enables us to determine the value of C2 for which a satisfactory 
solution exists. This value of C2 is interpreted as the chiral condensate. The numerical 
value of C2 itself is different for the D5 and the D7 solution. 

For the solutions we have found numerically we can then calculate the energy den- 
sities which (for the static configurations that we have been considering) are propor- 



tional to the negatives of the Routhians (19), (36) and (37). Finally, by determining 



the total energy of the D5 and D7 brane solutions corresponding to a given value of 
/ and V we can determine which configuration is the stable onej^ 

Numerical solutions for the profile of '0(r) for the D5 and D7 branes when u = 1 
and when / = 1 and / = 10 are plotted in figures [5] and |6] respectively. The case 
u = 1 is the one where the D7 brane is allowed to have a Minkowski embedding which 
is clearly seen in the numerical plot. For both of those cases, the D7 brane has lower 
energy than the D5 brane and it is therefore the preferred solution. 

We examine the solutions in the vicinity oi u = 1^ (u = 1.1,/ = 1) and {u = 
0.9, / = 1) in figures [T] and |8| respectively. In those cases, like the D5 brane, the D7 
brane must have a black hole embedding and this is indeed what is found numerically. 
The D7 brane has a long spike which goes to the Poincare horizon. Note from the 
plots that the D7 spike is significantly narrower than the D5 spike. The D7 brane 
then loses some of its energy advantage for small values of ln(r). It nevertheless still 
has lower total energy than the D5 brane, meaning that even for the non-gapped 
non-quantum Hall states, it is the preferred solution. 



The numerical solution for (z/ = 2, / = 10) is shown in figure 10, This is an integer 
filling fraction but the state is not a quantum Hall state. We shall explain shortly 
how to find the quantum Hall state with u = 2 (and which will be the energetically 
preferred one). In the meantime, we note that comparison of energies shows that this 
intermediate D7 brane state is preferred over the D5 brane. 



^When integrating the energy densities we use the earlier mentioned series expansions to esti- 
mate the contributions from the region of very large r and we approximate the densities with their 
asymptotic constant values to estimate the contribution from the small-r region. 
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Fi gure 5: Numerical solution with filling fraction v = 1 and / = 1. On the left-hand- graph, 
the profile of the function is plotted on the vertical axis versus ln(r) on the horizontal axis. The 
D7 brane is the dashed curve which comes to the horizontal axis. The D5 brane is the solid curve 
approaching the the horizontal axis at the Poincare horizon which is at the asymptotic left-hand- 
side of the graph. The second graph is the energy density on the vertical axis versus ln(r) on the 
horizontal axis. It should be integrated over ln(r) to get the total energy. Numerical computation 
shows that the D7 brane has lower energv, ^5-^7 = 1.71- 5 2 , and is therefore the preferred 

(27r) 2 Va 

state. 




-1.0 -0.5 0.5 1.0 -0.6 -0.4 -0.2 ^ 0.2 



Figure 6: Numerical solution with filling fraction u = 1 and / = 10. On the left-hand-graph, 
the profile of the function ip is plotted on the vertical axis versus ln(r) on the horizontal axis. The 
D7 brane is the dashed curve and the D5 brane is the solid curve. The second graph compares the 
energy densities. The D7 brane is the stable solution. Numerically, E5 — E7 = 2.34 • — — . 
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Figure 7: Numerical solution with filling fraction v = 1.1 and / = 1. On the left-hand- graph, the 
profile of the function ip is plotted on the vertical axis versus ln(r) on the horizontal axis. The second 
graph is the energy density on the vertical axis versus ln(r) in the horizontal axis. Like the D5 brane, 
in this case the D7 brane also approaches the Poincare horizon at the asymptotic left-hand-side of 
the graph. The D7 brane has lower energy. Numerically, E5 - E7 = 1.66 • B2 

(27r) 2 ^/A 




-12 -10 -8 -6 -4 -2 2 -3 -2 -1 



Figure 8: Numerical solution with filling fraction u = 0.9 and / = 1. On the left-hand-graph, 
the profile of the function ip is plotted on the vertical axis versus ln(r) on the horizontal axis. The 
second graph is the energy density on the vertical axis versus ln(r) in the horizontal axis. It should 
be integrated over ln(r) to get the total energy. The D7 brane has lower energy. Numerically, 
E5-E7= 1.43 • Bi. 

(27r)2 
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-3-2-1 1 -1.4 -1.2 -1.0 -0.8 -0.6 -0.4 -0.2 



Figure 9: Numerical solution with filling fraction v = 0.2 and / = 10. On the left-hand-graph, the 
profile of the function ip is plotted on the vertical axis versus ln(r) on the horizontal axis. The second 
graph is the energy density on the vertical axis versus ln(r) in the horizontal axis. In this case, the 
D5 brane has lower energy and is the stable solution. Numerically, E5 — E7 = —1.34 • — — 

(27r) 2 ^/x 




-8 -6 -4 -2 2 -1.5 -1.0 -0.5 



Figure 10: Numerical solution with filling fraction u = 2 and / = 10. On the left-hand-graph, 
the profile of the function ip is plotted on the vertical axis versus ln(r) on the horizontal axis. The 
second graph is the energy density on the vertical axis versus ln(r) in the horizontal axis. It should 
be integrated over ln(r) to get the total energy. The D7 brane has lower energy and is the stable 
solution. Numerically, E5 - E7 = 2.10 • . 
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Figure 11: Numerical solution with filling fraction u = 10 and / = 10. On the left-hand- graph, 
the profile of the function ijj is plotted on the vertical axis versus ln(r) on the horizontal axis. The 
second graph is the energy density on the vertical axis versus ln(r) in the horizontal axis. The D7 
brane has lower energy. Numerically, E5 — E7 = 0.22 • — — . 

(27r) 2 




Figure 12: Numerical solution which compares the energy densities of a (i^ = 2, / = 12) state, 
which has only a black hole embedding, with the energy density of two (i^ = 1, / = 6) states which 
have Minkowski embeddings. The vertical axis is energy density and the horizontal axis is ln(r). 
The energy density (per unit ln(r)) of the 2x (jy = 1,/ = 6) state is the dashed curve whereas 
the energy density of the (i^ = 2, / = 12) state is the solid curve. The sum of the energies of 
the two (i^ = 1,/ = 6) branes is lower that that of the {u = 2,/ = 12) brane. Numerically, 
^(n,/) = (l,12) - 2£;(^j) = (i,6) = 2.43 • .^JfvA^^* 
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Figure 13: This is a plot of (a constant minus) the total energies of two D7 branes, one of which 
has (i^ = 1, /i and the other of which has (i^ = 1, 12 — /i) for various values of /i between and 6. 
This clearly shows that the symmetric configuration where /i = //2 has the lowest energy and is 
the prefered state. 



We have done a numerical exploration of the region of smaller values of i^. We 
find that for > |, the D7 brane is stable and for < \ the D5 brane is the stable 
solution when / is in the region between one and ten. In figure[9]we plot the numerical 
solution for {v = 0.2,/ = 10). For that solution, we find that the D5 brane is the 
configuration with lower energy. 



Some results for a larger value of u are shown in figure [10] and [TT] where we 
display the solutions with (z/ = 2,/ = 10) and {u = 10, / = 10). In both of these 
cases, the D7 brane has lower energy. We have also looked at an extreme case where 
{u = 100, / = 100) and we find that, even though the energies become very close, D7 
is still preferred. 

Then, we consider some positive integer filling fractions and compare the following 
solutions, (z/ = 2, / = 12) and two D7 branes, each of which have (z/ = 1, / = 6). The 
energy densities of the two cases are plotted in figure [9l The energy of the two branes 
with z/ = 1 is lower than the single brane with u = 2. This can be attributed to the 
fact that the former has a Minkowski embedding whereas the latter has a black hole 
embedding. The Minkowski embedding has an energy advantage in the small ln(r) 
region. 

We also consider asymmetric splittings of the A^5 D5 branes into two D7 branes 
when u = 2. For / = 12 we consider a number of cases where the first D7 has u = 1 
and fi and the second D7 has u = 1 and /2 so that /i + /2 = /• The energies for a 
number of cases with /i varying from 2 to | = 6 are plotted in figure [isj 

Finally, we examine the question of whether the chiral symmetry restoring phase 
transition which occurs for the D5 brane as the filling fraction u is increased to its 
critical value Uc ^ 1-68/ given in equation (20) is seen when the D5 brane is blown 
up to a D7 brane. In the above discussion, we have learned that in a certain region 
of the (z/, /) plane, the D7 brane has lower energy than the D5 brane and is therefore 
the preferred solution. Now we will report preliminary investigation of the boundaries 
of the region where it is more stable. First of all, we have seen that, for moderate 
values of / (in the range 1-10), the D7 brane takes over from the D5 brane when u 
becomes larger than ^, at least in the region where the D5 brane would have a chiral 
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Figure 14: The solution for a D7 brane where u = 1 and / = 0.1 is depicted (dotted hnes) in the 
left-hand-graph and its energy density (as a function of ln(r)) in the right-hand graph. The solid 
line in both graphs is the profile and energy density of the constant tjj = ^ solution which is the 
only D5 solution in the u > Uc regime. For this value of /, Uc ^ .168 and u = 1. By integrating 
the energy density, we clearly see that the chiral symmetry breaking D7 brane is still energetically 



Figure 15: The solution for a D7 brane where u = 0.9 and / = 0.1 is depicted (dotted hnes) 
in the left-hand-graph and its energy density (as a function of ln(r)) in the right-hand graph. The 
solid line in both graphs is the profile and energy density of the constant = f solution which is 
the only D5 solution in the v > Vc regime. For this value of /, ^ -168 and v = 0.9. By integrating 
the energy density, we clearly see that the chiral symmetry breaking D7 brane is still energetically 



preferred, — = 1.11 • 





preferred, E^ — Ej = 1.26 • 



(27r)2 
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symmetry breaking solution with non-constant ip{r). Thus, when < |, at least in 
the regime that we can analyze, the standard D5 brane phase diagram applies. 

The D5 embedding equation depends only on the ratio j. If we consider the 
conventional picture of the D5 brane with a magnetic field and charge density, at 
very small j, it is in a chiral symmetry breaking state. The function ip^r) has a 
non-constant profile which goes to | as r ^ oo (as in equation (22) ) and goes to zero 
as r ^ 0. We are considering only the case where Ci = and then, for a non-constant 
'0(r), C2 must be nonzero. If we increase the ratio j the magnitude of the chiral 

condensate C2 decreases and becomes zero when Ty/f = At this critical point, 

there is a phase transition with Berezinski-Kosterlitz-Thouless scaling [12]. On the 
other side of the critical point, where / f > the only solution of the D5-brane 
theory is the constant = ^. 

First, we observe that the D7 brane equation also has the same constant solution 
as the D5 brane equation = | with the same energy and it becomes unstable at 

the same value of j. This means that, whenever v j f > the constant = | is at 
least a meta-stable solution for the D7 brane. 

However, the D7 brane solutions with non-constant '0(r) have a different behavior 
than those for the D5 brane. First of all, they depend on both parameters v and 
/, rather than just the ratio j. Secondly, numerical computations show that they 

continue to exist, even in the region j > ^ 1.68 where the constant solution 
= I is a competitor. An example with = 1 and / = 0.1, so that y = 10, is shown 



in figure [14} Not only does the solution exist there, but it has lower energy than 
the constant solution and it is therefore the energetically preferred configuration. 
Ungapped chiral symmetry breaking states are also energetically favorable in this 



region. Figure 15 displays an example where j = 9 and the D7 brane has a black hole 
embedding anoTt is also stable. In an attempt to find solutions which restore chiral 
symmetry, we have studied some examples with extreme values of j where the D7 
brane still has lower energy. Indications are that, for a very large range of parameters 
in the region > |, the chiral symmetry breaking D7 brane is the lowest energy 
solution. In the above, we argued that the constant solution should be metastable 
in the whole region > ^ 1.68. The existence of non-constant solutions with 
lower energy there implies that the constant and non-constant solutions are separated 
by an energy barrier. This in turn implies that, if their energies became equal or if 
the constant solution became energetically favorable in some region of the parameter 
space, the chiral symmetry restoration phase transition would be first order. At this 
point, we have not found numerical evidence for such a phase transition, however, 
this point certainly deserves additional study. 

Recall that our Hall states for = 1, 2, . . . , A^5 have D7 branes with /, //2, f /N^ 
(assuming that the A^5 D5 branes are shared symmetrically between the D7 branes) 
and therefore have j = 1, 2, . . . , A^5. There is still the open question as to how large 
this ratio can be before chiral symmetry is restored, or alternatively, how many of 
these quantum Hall states are stable. As we have discussed above, we have found a 
numerical solution where j = 10 (in figure flil) where the chiral symmetry breaking 
Hall state is still the one which is energetically prefered. We have also examined the 
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extreme case of u = 1 and / = 0.01, that is ^ = 100 and it was also stable. 



5 Conclusion 

We have formulated a mechanism by which integer quantum Hall states can appear 
in AdS/CFT holography which uses probe D5 branes to construct strong coupling 
A/" = 4 Yang-Mills theory with fundamental representation matter fields occupying 
a supersymmetric 2+1-dimensional defect. The external magnetic field and finite 
charge density that are needed in order to obtain a Hall state break supersymmetry 
(even in the free field limit, the fermions and bosons have different Landau level 
spectra). We propose a number A^5 of D5 branes blown up to form a D7 brane as a 
candidate for the stable solution over a parameter range which we can easily analyze 
when the number of D5 branes is large, A^5 ^ V^. 

We find that for a range of filling fractions u > 0.5 and for / = ;^A^5 of order 
one-to-ten, the D7 brane has lower energy than the D5 brane and it is therefore the 
preferred solution. When u < 0.5 the D5 brane becomes stable, at least for the range 
of / that we consider. Some of the solutions with their energy densities are shown in 
figures [5p2 



When u > 0.5 we find stable D7 brane solutions with non-constant '0(r) for a large 
range of values of the parameter /, including examples in the region where j > 
This value of y is a critical point of the D5 brane where there is a phase transition 
which restores chiral symmetry. The situation is very different when u > 0.5 where 
the D7 brane is more stable. In that case, the chiral symmetry breaking persists to 
larger values of j. We have not been able to determine whether there is a critical value 
of this parameter. The competing chirally symmetric solution seems met ast able in 
this region, indicating existence of an energy barrier between the solutions. Therefore 
one would expect a chiral symmetry restoring phase transition to be of first order, if 
it exists at all. 

We have used a comparison of energies to decide which D brane configurations are 
more stable. However, we have not done a full stability analysis. Particularly at finite 
density, instabilities to formation of non-translation invariant phases such as striped 
phases do occur for the D7' system [46j and such instabilities could also appear in the 
present system. We have left checking this important possibility for future work. 

The mechanism by which we get higher filling fraction Hall states is by forming 
multiple D7 branes which share the total charge density and which each have unit 
filling fraction u = 1. In principle, states which have filling fractions = 1, 2, . . . , A^5 
could exist. How many of them exist is a dynamical question. We have found nu- 
merical solutions where u is as large as 10 and the Hall state still has the lowest 
energy. 

We have used numerical computation to show that, for the first higher filling 
fraction state = 2, the two D7 branes with u = 1 are stable in that they have lower 
energies than the single D7 brane with = 2, which in turn has lower energy than 
the D5 branes with these parameters. The main reason for their stability is that they 
are Minkowski embeddings, so they gain energy by not having spikes going to the 
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Poincare horizon. Indications are that, in the multi D7 brane solution, at least for 
z/ = 2, the state with D5 branes evenly distributed between the two D7 branes is 
favored over asymmetric distributions. The v D7 branes are therefore coincident and 
have an SU(i^) symmetry. However, more work on mapping the energy landscape of 
multiple D7 branes is definitely needed before this is a firm conclusion. 
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